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Abstract 

A (n + l)-dimensional cosmological model with a set of scalar fields 
and antisymmetric (p+2)-form is considered. Some of scalar fields may 
have negative kinetic terms, i.e. they may describe "phantom" fields. 
For certain odd dimensions (D = Am + 1 = 5, 9, 13, ...) and (p + 2)- 
forms (p = 2m — 1 = 1,3,5,...) and non-exceptional dilatonic coupling 
vector A we obtain cosmological- type solutions to the field equations. 
These solutions are characterized by self-dual or anti-self-dual charge 
density forms Q (of rank 2m) and may describe the maximal set of 
branes (i.e. when all the branes have non-zero charge densities). Some 
properties of these solutions are considered, e.g. Kasner-like behavior, 
the existence of non-singular (e.g. bouncing) solutions and those with 
acceleration. The solutions with bouncing and acceleration take place 
when at least there is one "phantom" field in the model. 
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1 Introduction 



In this paper we investigate composite electric S'-brane solutions (space-like 
analogues of D-branes) in an arbitrary number of dimensions D with a set 
of scalar fields and (p + 2)-form. The (p + 2)-form is considered using a 
composite electric ansatz and the metric is taken as diagonal. All ansatz 
functions for the metric, form field and scalar field are taken to be dependent 
upon one distinguished coordinate u. (For S'-brane solutions see [Tj- [TO] and 
references therein.) In [I]g\ a special case of this model containing one scalar 
field with positive-definite kinetic term was considered. 

Here, as in [16], we investigate our system similar to the approach used 
in ^7j. We use the same sigma- model approach with constraints and use 
the result of the previous work [T^], where it was shown that it is possible 
to satisfy the constraints maximally {i.e. all the branes may carry non-zero 
charge densities) in certain odd dimensions D = 4m + 1 = 5, 9, 13, .... 

Here we also examine cosmological-type solutions in these odd dimen- 
sional cases, and discuss some of their interesting features such as their 
Kasner-like behavior, existence of non-singular (e.g. bouncing) solutions and 
solutions with acceleration among them. 

In the next two sections we will give the set up for the system of D = 
n + 1 dimensional gravity with a set of scalar fields and (p + 2)-form field. 
For the conditions considered in this paper (diagonal metric, composite Sp- 
brane ansatz for the antisymmetric (p + 2)-form field, and all the ansatz 
functions depending only on one coordinate) this system can be reduced to 
a 1-dimensional a-model with contstraints. 

In section 4 we consider the quadratic constraints for the charge densities 
of the branes. Here we briefly remind one of the main result of where it 
was found that these constraints have "maximal" solutions with all non-zero 
brane charge densities in particular odd dimensions with particular ranks of 
form fields: D = 4m + 1 and p = 2m — 1, respectively (m = 1,2,...). 

In section 5 we investigate cosmological solutions to the field equations for 
these odd dimensions. We look at the proper time behavior of the simplest of 
these solutions. We show that certain solutions exhibit Kasner-like behavior 
at the early times, and some them have also such behavior for big times. 
Also we find an accelerated expansion of n-dimensional space for certain 
negative values of dilatonic vector coupling squared, i.e. when there at least 
one "phantom" scalar field in the model. 

In section 6 we single out two subclasses of non-singular solutions, e.g. 
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with bouncing behavior of the scale factors. Here we also show that solutions 
with bouncing take place only when there is at least one "phantom" scalar 
field in the model. 



2 D-dimensional gravity coupled to g-form and 
scalar field 



Here we consider the model governed by the action 



S= I d D z^/\g~\ 

' M 



R[g] - h ab g MN d M ^ a d N ^ b - ^exp(2A a ^)F 2 



(2.1) 



where g = gMNdz M ® dz N is the metric, (p = (ip a ) is a set (vector) of scalar 
fields, a = 1, . . . , I, A = (A a ) G R is a constant vector (set) of dilatonic 
couplings, (h a },) is a symmetric non-degenerate I x I matrix and 

F = dA= -F Ml Mq dz Ml A ... A dz M \ (2.2) 
ql 

is a g-form, q = p + 2 > 2, on a D-dimensional manifold M. Without loss of 
generality the matrix (h a b) may be considered as a diagonal one. 
In (12. 1|) we denote \g\ = \ det(gM7v)l an d 

F 2 = F Ml ... Mq F Nl ... Nq g M ^ . . . g M « N «. (2.3) 

The equations of motion corresponding to (j2.1J) are following 

Rmn — ^9mnR = Tmn, (2.4) 
A[^ a -^e 2 V = 0, (2.5) 

V Ml [g}(e 2 ^F Ml - Mq ) = 0, (2.6) 

where A a = h ab Xb, (h ab ) is a matrix inverse to (h a b) and \(p = \ a tp a - 

In fl2.5j) and ()2.6|) . A [g] and X/[g] are Laplace-Beltrami and covariant 
derivative operators corresponding to g. Equations ()2.4|) . (j2.5j) and ()2.6|) are, 
respectively, the multidimensional Einstein-Hilbert equations, the "Klein- 
Fock-Gordon" equation for the scalar field and the "Maxwell" equations for 
the g-form. 
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The stress-energy tensor in (|2.4|) can be split up as 

Tmn = 

with 

T MN [<f,g) = h ab (d M ^ a d N ip b - -g M Nd P (p a d p (p b ), 



Tmn[F,9] = -r 



--QmnF + qFuM-2. ...Mq-^N 



M 2 ...M q 



(2.7) 

(2.8) 
(2.9) 



being the stress-energy tensor of scalar fields and g-form, respectively. 
Let us consider the manifold 

M = (u_,u+) x R n 

with the diagonal metric 

n 

g = we 2 ^ u) du <g> du + eW{u )e i d V* ® d v\ 

i=l 



(2.10) 



(2.11) 



where w — ±1, and u is a distinguished coordinate. The functions ^(u), 4> % {u), 
the scalar fields <^ a (w) and the g-forms are assumed to depend only on u. 
Here 

Si = ±1 (2.12) 

are signature parameters, i = 1, . . . ,n. When w = —1 and all Si = 1 the 
solutions are cosmological. The functions 7,0*: («_,«+) — > M are smooth. 
Here as in JH] we define 



fi o = {0,{l},...,M,{l,2},...,{l,...,n}} 
which is the set of all subsets of 

I Q = {1, . . . ,n}. 



(2.13) 



(2.14) 



For any I — {ii, . . . , ik} € Qq with i\ < . . . < i^, we define a form of rank 
d{I) = k 

t(I) = dy h A ... A dy ik , (2.15) 

The corresponding brane submanifold is described by coordinates y 
We also define the ^-symbol as 



^1 oi^k 



(2.16) 
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Here we adopt the following electric composite Sp-brane ansatz for the 
field of the (p + 2) -form 



F = ^ d& A r(J) 



(2.17) 



where the set 



Q e = {I e n \d(I) = q-l=p+l} 



(2.18) 



contains all subsets of flo of the "length" p + l, i.e. of the form ii, i p }. 

We assume that the scalar potential and all scalar fields depend on the 
distinguished coordinate 



3 cr-model representation with constraints 
3.1 cr-model 

As it was shown in J7| (see Proposition 2 in [TJj) the diagonal part of Einstein 
equations ()2.4j) and the equations of motion ()2.5|) - ([2.6|h for the ansatz given 
in (|2.1ip . (|2.17|) - ()2.19|) . are equivalent to the equations of motion for a 1- 
dimensional cr-model with the action (see also [HI El) 



(2.19) 




(3.1) 



the overdots represent differentiation with respect to u, i.e. 4-. 
The factor M is the lapse function given by 



M = exp(7o — 7) > 



(3.2) 



where 
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(3.3) 



i=l 



The factor in the exponent is given by 




(3.4) 
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Here and it what follows 

Gij — Sij — 1 (3-5) 

are components of the "pure cosmological" minisupermetric matrix, i,j = 
l,...,n, see [22 123- 

The generalized "Maxwell equations" of ()2.6j) 

I {°M-1U'W) = (3.6) 

can be readily integrated to give 

$ 7 = Q(/)exp(2f/ / ), (3.7) 

where Q(I) are constant charge densities, / G Q e . 
We will fix the time gauge to be harmonic one 

7 = 70, M=l. (3.8) 

We also introduce collective variables x = (x A ) = ((f) 1 , ip a ) and a "trun- 
cated" target space metric 

G = G A Esdx A <g> dx B = Gijdcj) 1 <g> d<jt + h ab dip a <g> dip b , (3.9) 

(C - ) = ( G o' J I)' (3 ' 10) 

The U 1 - vectors defined in ()3.4|) can be written as U 1 \x) = U I A x A with com- 
ponents 

(U I A ) = (5 iI ,-X a ) , (3.11) 

where 

is an indicator of i belonging to J. For fixed charge densities Q(I), I <E Q e , the 
equations of motion for the cr-model in (|3.1|) are equivalent to the Lagrange 
equations corresponding to the Lagrangian 

Lq = \g ab x a x b - V Q , (3.13) 

with the zero-energy constraint 

Eq = -G AB x A x B + V Q = (3.14) 



imposed. Here 

V Q = \Y, S ^Q 2 W exp[2t/ 7 (x)]. (3.15) 

In section 5 we will find explicit solutions to the field equations resulting 
from (j3.13j) - (|3.15|) for certain dimensions D = 4m + 1 and positive 

3.2 Constraints 

Due to diagonality of the Ricci-tensor for the metric (|2.11|) the non-diagonal 
part of the Einstein equations ()2.4j) reads as follows 

!fy = 0, i^j. (3.16) 

This leads to constraints on the charge densities Q(I). The non-diagonal 
components of stress-energy tensor are proportional to 

e 2 ^F lM2 ... Mq F 3 M *- M \ (3.17) 

with i j. From ()2.17|) ()3.4j) and ()3.7|) we obtain for the (p + 2)-form 

F = . —r: Qi ...i p exp(20 4 ° + . . . + 2<P* - 2\0)du A dy<° A ... A dy ip (3.18) 
(p + ly. 

Inserting this in (J3.17|) we are led to the following constraint equations on 
charge densities [TT] 

n 

^ij = ^ ^ Qii\...i p Q • ■ • 0; (3.19) 

il,...,i p =l 

where i ^ j; i,j — 1, . . . , n. (Tjj is proportional to exp(— 2X0 — 27 + 2<p % + 
2o l )C, l lor//./.) 

Here p = q — 2 and Qi i 1 ..i p are components of the antisymmetric form of 
rank p + 1 = q — 1 and 

Qi h...i p = Q({io, h, ■ ■ ■ , i P }) (3.20) 

for i < ii < . . . < i p and {i , i 1 , . . . ,i p } e Q e . 

The number of constraints in ()3.19|) is n(n — l)/2. It was shown in JH] 
that these constraints can be satisfied when the dimension of space-time takes 
certain odd values (see next section). 
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4 Solution to constraints for D= 4m+l 

As in jTHj we rewrite the constraints ()3.19|) as 

n 

h,...,i p =l 

i 7^ j'i h J ' = 1j • • • 5 n ! where we introduce "running" charge densities 

p 

0i «...%, = Qia*i...ip JJexp(0 ife ). (4.2) 

fc=0 

Here the indices are lifted by the flat metric 

77 = sidy 1 ®dy l + ...+ e n dy n ® rfy n = 77 a6 tft/ a ® dy b , (4.3) 

where (r? afc ) = (r/ ab ) = diag(ei, . . . , e n ), i.e. Q l °--' l p = 77*0*0 . . ■r] lpl ' p Qi' ...i' p - Here 
we denote Cf = Cf exp(0* + ft) with C| = C ik r] kj . 

It was shown in JE] that when the total space-time dimension is the 
following one 

D = n + 1 = 2(p + 1) + 1 = 4m + 1 = 5, 9, 13, . . . (4.4) 

with 

p = 2m- 1 = 1,3,5,... (4.5) 

and 

e 1 ...e n = l (4.6) 
then any self-dual or anti-self-dual "running" charge density form 

Qi ...i p — i J- — ^y £io...ipj ...j p Q :>0 ]p — i(*Q)i — ip m (4-7) 

satisfies the constraint equations f)4.1j) (or, equivalently, eqs. (|3.19j) ). Here 
the symbol * = *[rj\ is the Hodge operator with respect to 77. The relation 
(I4.fjjl means that the number of time-like coordinates is even. 

Remark. It was shown in ^Bj that in the special case: D = 5, p = 1 
and £1 = . . . = £4 = 1 the solutions presented above are the only solutions 
to constraints ()4.1j) when all (* < j) are imposed to be non-zero ones. 
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5 Cosmological-type solutions for D= 4m+l 

Here we consider cosmological type solutions for the dimensions from (|4.4|) 
when all ^-coordinates are space-like, i.e. 

e 1 = ... = e n = l. (5.1) 

It was shown in jTH] that the definition of running constants (|4.2|) and 
the (anti-) self-duality of the charge density form (J4.7|) imply that all scale 
factors are the same up to constants: 

ftfa) = (f)(u) + c\ (5.2) 

Without loss of generality we put c 1 = which may always be done via a 
proper rescaling of y-coordinates. This also implies that non-running charge 
density form Qi ...i p is self-dual or anti- self- dual in a flat Euclidean space M. n , 
i.e. 



Qi Q ...i p — i , i\, £ io---ipjo-- jpQ :>0 3P — ^{*Q)io...i p - (5-3) 



1 

(j, . 1 ""••'<••' •'<" 

The Lagrangian and total energy constraint are given by 



L Q = ^-^q + ^WV, (5-4) 

E Q = ^#' + ^ + ^WV = 0, (5.5) 
with the potential being 

V Q = \y (^f/)(-xpl 2> <> k -2\?\. (0.6) 



^Q 2 (/)exp (2^/ -2A^) 
/ V fee/ / 



see dsnsj-dsnnj) 

The field equations for and ip a from Lq are the following ones 

^G^ + ^Q 2 (/)5}exp(2^0 fc -2A^) =0, (5.7) 
3=1 i \ kei / 

/ V fee/ / 
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a = 1, . . . , I. 

Since all 4> % = <p we get YTj=x = Z)™=i(^i _ l W = 0- ~ 

As in |16j . denning 

Q 2^J2Q 2 (I)^0 (5.9) 

and noting that due to (|5.3I) 

for any i = 1, n, we could write the field equations (j5.7J) and ()5.8|) in the 
following form 

' Q 2 exp(n(f)-2\0), (5.11) 



2(n-l, 

y3 a = A a Q 2 exp(n0-2A^), (5.12) 

a = 1, . . . ,1. 

Equations (|5.11jl and ()5.12|) imply 

f = -2Ae 2f , (5.13) 

where the definitions 

f = -<j>-\& (5.14) 

A = ^-K, K = \ 2 --^- (5.15) 
2 4(n — 1) 

are assumed. Here and it what follows A 2 = h a b\ a X b . 
The first integral of (J5.13)) is 

1 -f 2 + Ae*f= 1 -C, (5.16) 

where C is an integration constant. 
Let K ^ 0, or 

Equation ()5.13|) has several solutions [16J: 

/ = -In [^|2v4| 1/2 ] (5.18) 
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with 



1 



C 
1 



sinh 



(u -«o)vC 



sin 



cosh 



U — Uq 
1 



C 



A < 0, C > 0; 

A < 0, C < 0; 
A < 0, C = 0; 
A > 0, C > 0. 



(it -iio)vC 

Using 1)5.11)) and (|5.12l we get the following relationship 

(p a = 2(n - l)A a 0, 

which has the solution 

<p a = 2(n - l)A a + C> + C£, 



(5.19) 

(5.20) 
(5.21) 
(5.22) 



(5.23) 



(5.24) 



where C% , are integration constants, a = 1, . . . ,1. Combining ()5.24)) with 

1 [/(«) + A(d 2 « + Ci) 



and 



2(1 -n)AT 
A" 



A" 



/(tt) + A(C r 2 « + C r i) + C> + C 1 a , 



(5.25) 
(5.26) 



0=1,...,/. 

Substituting this to the zero-energy constraint 

E Q = in(l - n)0 2 + -/i^V + ^Q 2 e 2/(n) = 
and using (J5.16)) we get 



or, equivalently, 

C = (AC 2 ) 2 - K(C 2 ) 2 = (\C 2 ) 2 - A 2 ((7 2 ) 2 + 



(5.27) 



(5.28) 



n 



4(n-l) 



(C 2 ) 2 . (5.29) 



where (C 2 ) 2 = /i a6 C 2 a C 2 ' and AC 2 = A Q C 2 a . 
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Thus, the solutions for the metric, scalar fields and (p + 2)-form read: 

n 

ds 2 = we 2n ^du 2 + e 2 ^^(<V) 2 (5.30) 

i=i 

^ = ~ [f(u) + X(C 2 u + C 1 )] +C«u + Cl (5.31) 
F = e 2f{u) duAQ, Q = t^~ 



j—yQi ...i p dy l ° A.-A^, (5.32) 



a = 1, . . . , Z, with given by (|5.25jl and the function f(u) given by (|5.18|) 
and dSHSl, IjOTjl . JE2U), ((£22]). Here the charge density form Q of rank 
n/2 = 2m is self-dual or ant i- self- dual in a flat Euclidean space M. n : Q = 
± * Q, the parameters Cf , C obey ()5.29|) and the vector of dilatonic coupling 
constants A is non-exceptional one, see (|5.17|) . 

5.1 Special attractor solution for = 

Here we consider the simplest cosmological type solution with Cf = and, 
hence, C — 0, i.e. the solution corresponding to ()5.21|) . For this solution 
A < and hence 

A 2 < A 2 . (5.33) 

Without loss of generality we put u = 0. 

Let us introduce a "proper time" coordinate: 

dr = -e n ^ u) du, (5.34) 

where 

1 ln^^Aj 1 / 2 ), A = Ae- 2XaC K (5.35) 



r 2(n-l)K 

Integrating ()5.34|) and taking a suitable choice of reference point we get for 
A 2 + -A 2 

\a\\2A\ x l 2 r = (u\2A\ 1/2 ) a , (5.36) 



where u > and 



a = £±*a. (5.37) 
A 2 - A 2 V ; 



We note that for A 2 > —A 2 , we get a < 0, and hence r = t{u) is mono- 
tonically decreasing from infinity (when u — 0) to zero (when u = oo). 



12 



The metric (|5.3U|) now reads 



ds 2 = wdr 2 + Bt 2v ^(dy*) 2 , (5.38) 

i=l 

where r > and 

v= * -, B = (\a\\2A\^r. (5.39) 

n + 4X z {n — 1) 

For scalar field we get the following relations 

<p a = 2\ a (n - l>lnr + pg, (5.40) 



(dr A Q) \2A\ 


(a-l)/2a 


2a\A\ 


T 


(la 


| r )l/a 



a = 1,...,/; where 5 = fi exp(-4A a ^g) and B = (\a\ \2A\^ 2 ) 2,/ . 
For (p + 2)-form we get 

a /T\ lo /I I fa— l)/2a 

exp(-2A ^g). (5.41) 

By putting A = and w = — 1 in the above solution we get a cosmo logical 
power-law expansion with a power u — 2/n that is the same as in the case of 
D = 1 + n dust matter with a zero pressure (see, for example [201 |2U [22J ) . 
This coincidence was explained in [IB], i-e. it was shown that the collection 
of branes with charge densities obeying (anti)-self-duality condition ()5.3j) 
behaves as a dust matter. 

The solution given by relations ()5.38jh (|5.38|) . (|5.41|) is an attractor solu- 
tion in the limit r — > +oo, or«-^ +0, for the solutions with A < given by 
(|5.19|) . This follows just from the relation sinhw ~ u for small u. 

Accelerated expansion. When 

v > 1 (5.42) 

we get an accelerated expansion of n-dimensional space. This relation is 
equivalent to the following one (n > 2) 

n < A 2 < - <?- 2 l (5.43) 



4(n-l) 4(n 

This means that A 2 < and, hence, the matrix (h a b) is not a positive-definite 
one. So, there should be at least one scalar field with negative kinetic term 
(i.e. so-called phantom field). 
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5.2 Kasner-like behavior for r — ► +0 

Now, we consider u — > +00 asymptotical behavior of solutions with i) sinh- 
and ii) cosh- functions corresponding to (j5.19|) and ()5.22j) . respectively. In 
both cases we find that the asymptotic behavior is Kasner like. We remind 
that in general case such behaviour is described by the following asymptotical 
relations for the metric and scalar fields [T81 ITU] 

n 

ds 2 as = wdr 2 + r 2al A t W) 2 , ¥4 = o£ In r + ^, (5.44) 
i=i 

where A4 > 0, v?o are constants. The Kasner parameters obey 

n n 

£V = ]T(«*) 2 + h ab a a v a b v = 1. (5.45) 
1=1 1=1 

In our case all a 1 are coinciding since we have an isotropic expansion 
described by the metric. 

Here we consider for simplicity the case when the matrix (/i a b) is positive 
definite one. In what follows we consider solutions with C > 0. We also put 
C 2 ^ that imply (C 2 ) 2 > 0. 

In the first case i) with sinh-dependence and K < (or, equivalently, 
A 2 < Aq) the proper time r is decreasing when u — > +00. From eqs. ()5.18|) . 
(|5.19|) , (|5.22|) , (|5.25|) and (|5.34|) we find the following asymptotic behavior as 
u — > +00 

~ —bu + const, b = — — — > 0, (5.46) 

2{n — 1)K 

v a ~ [C 2 a - 2(n - l)b\ a }u + const, (5.47) 
t ~ const exp(— n6ti). (5.48) 

The inequality (|5.46|) follows just from ()5.29j) . 

Using these asymptotic relations and writing everything in terms of proper 
time one find that the metric and scalar field take the following asymptotic 
forms 

n 

ds 2 as = wdr 2 + T 2 / n A Y,W)\ (5.49) 

i=i 

^ s = ^lnr + ^, (5.50) 
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as r — > +0. Here A > 0, y?Q are constants and 

2(n - l)A a C% 



n n6 
It may be verified that the Kasner relation 



< = i " rf- ( 5 - 51 ) 



= 1 - - (5.52) 



1 



71 



is satisfied identically (in our case all a 1 = 1/n). 

In the second case ii) with cosh-dependence and K > (or, equivalently, 
A 2 > Ag) the proper time r decreases as u — > +oo for A a C*2 > and increases 
for A a C*2 < 0. In this case we also get an asymptotical Kasner type relations 
lj05)) - (pr5Dl) boths in the limits r -> +0 and r -> +oo. 

The Kasner parameters for scalar fields have the following form 

= Iz-W _ CI 
n no T 

where 

T 2(n-l)K v ' 

The signs "-" and "+" refer to the asymptotical behaviors for r — > +0 and 
r — > +oo, respectively. 

In both sinh — and cosh — cases the Kasner set a = (a 1 = 1/n, a" ) 
corresponding to the limit r — > +0 obeys the inequalities 

U I (a) = -X a a a v + J2 al = - X -K + l >0 ( 5 - 55 ) 

16/ 

for all brane sets /. This is in agreement with a general prescription of the 
billiard representation from [18J. 

In the case ii) we obtain that the Kasner set corresponding to the limit 
r — > +oo obeys the inequalities 

U\a) < (5.56) 

for all /. 

Relations (|5.55J) and ()5.56|) are nothing more then the conditions for van- 
ishing of the g-form term exp(2\ a (p a )F 2 in the Lagrangian (j2.1jl in Kasner 
regimes with r — > +0 and r — > +oo, respectively (see [XHJ ) - 
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Scattering law for Kasner parameters. The Kasner parameters for 
in" (r — > +0) and "out" (r — > +00) asymptotics are related as follows 

2(n-l)X a a a - 2(n ~ 1)Aa 

— £M ^ „ (5>5?) 

This relation could be readily obtained from formulas (|5.5Hjl and (|5.54|) . 
Let us explain the geometrical sense of (|5.57|) putting for simplicity h a b = 5 a b- 
Vectors o.^- and belong to the (/ — l)-dimensional sphere of radius 

R = — i. Let us introduce a vector A = H^zl)_ The endpoint of this 
vector is out of the sphere (of radius R) due to K > 0. If we consider a line 
connecting the endpoints of the vectors A and «,,_, then we get another point 
of intersection of this line with the sphere. This another point corresponds 
to the vector a Vt+ . 

Remark. In we have suggested an example of never-ending oscil- 
lating behavior near the singularity supported by scalar fields with multiple 
exponential potential. In ^J] we have found a "reflection law" for scalar 
Kasner parameters using the general relation from Here we suggest an 
example of "Kasner scattering", when Kasner behaviour takes place both 
"in" (r — > +0) and "out" (r — > +00) asymptotics. In this case the "scatter- 
ing" also takes place only for scalar part of Kasner parameters, while Kasner 
parameters for the metric are coinciding in both asymptotics. 

Remark. Here as in we also deal with the illumination of a sphere 
by a point-like sources. Indeed, if we put a source of light into the endpoint 
of the vector A then we get two domains on the sphere: i) the domain 
illuminated (strongly) by the source of light; and ii) the shadow domain. 
Due to inequalities (j5.55J) and ()5.56|) the vector a^- belongs to the shadow 
domain and the vector 5^+ belongs to the illuminated one. The map (|5.57J) 
transforms from the shadow domain into another one. 



6 Non-singular solutions 



Here we single out subclasses of non-singular solutions, e.g. with bouncing 
behavior of scale factor 



exp 



1 

2(1 -n)K 



(6.1) 
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where 

/ = / + A«(C£u + C?). (6.2) 

In this section we put 

A 2 > -A 2 . (6.3) 

Let us consider the case K < 0, or equivalently, A 2 < A 2 ,- In what follows 
we use synchronous or proper time variable from (J5.34j) . or, equivalently, 

du'e n ^ u '\ (6.4) 

For C < we have a bouncing behavior of scale factor 
It follows from (j6.1J) . ()6.2|) and ()6.4|) that the solution is a non-singular 
one for r G (— oo, +oo) in such three cases 

(A+) A a C 2 a > VC, C>0; (6.5) 
(4>) \ a C% = Vc, C > 0; (6.6) 

C < 0. (6.7) 

Indeed, for u = uo, the integral in ()6.4j) is divergent in all three cases due 

to 

e n ^ u) ~ ( u - u ) a -\ (6.8) 

as u — > uo + 0, since a < (see, (j5.37}l ) according to ()6.3j) . This implies 
r(«) — > +CX3 as u — > uq -\- 0. The relation ()6.8|) just follows from relations 
sinh(x) ~ sin(x) ~ x for small x. 

In the case (^4_) the function <p{u) is defined on the interval (uo, u + ) and 

e "*(«) ^ ( u+ - M )«-\ (6.9) 

as it — >• n + — 0. This follows just from sin(7r — x) = sin(x) ~ x for small x. 
Relation ()6.9j) implies r{u) — > — oo as u — > u+ — 0. Thus, in the case (A-) 
the solution is defined for r G (— oo, +oo). 

For (A + )— and (A Q )— cases the solution is defined for u G (u Q , +oo). 
When C > the function (|6.2j) has the following asymptotical behavior 

f (u-u )VC + \ a C%u + const (6.10) 

for u — > +oo, that implies the divergency of the integral in ()6.4)1 for u = +oo. 
An analogous consideration could be done for C = (with logarithmic term 
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appearing instead of the first term in r.h.s. of (|6.1U|) ). Thus, in both cases 
(A + ) and (A ) the solution is defined for r G (—00, +00). 

In the cases (A + ) and {A J) we also get a bouncing behaviour of the 
scale factor a(r) at some point Tb, i.e. the function a(r) is monotonically de- 
creasing in the interval (—00, r^) and monotonically increasing in the interval 
(r 6 , +00). 

In the case (A ) the solution is non-singular, the function a(r) is mono- 
tonically increasing to infinity from a non-zero value. 

Remark. A similar behavior takes place in D-dimensional model with 
two sets of scalar fields: "usual" (p and phantom ip ones and exponential 
potential depending upon ip, see |23] . 

Let us show that the bouncing behavior for K < takes place only when 
the matrix (h a b) is not positive definite, i.e. when at least one scalar field is 
a "phantom" one. Indeed, in the case (A + ) we get squaring ()6.5|) 

(X a C a 2 f >C= (\C 2 ) 2 - K(C 2 ) 2 , (6.11) 

and, consequently, (C 2 ) 2 = h ah C 2 C\ < 0. The same inequality we get in the 
case (AJ) just from the relation C = (XC2) 2 — K(C2) 2 < 0. This means 
that the matrix (h a b) is not positive-definite, i.e. we have a certain number 
(k > 1) of "phantom" fields among (p a . 

Let us consider the case K > 0, or equivalently, A 2 > A 2 ,. We get from 
(fn~221) C > 0. 

The solution is a non-singular one for r G (—00, +00) in such two cases 

(B + ) \\ a CI\<VC; (6.12) 
(B ) \\ a C«\ = VC. (6.13) 

Indeed, for K > the solution is defined for u G (—00, +00) and the 
function ()6.2j) has the following asymptotical behavior 

/ ~ T ( U _„ )VC + A tt C> + const (6.14) 

for u — > ±00. This relation implies the divergency of the integral ()6.4j) for 
u = ±00, when restrictions (-B+) and (B ) are imposed. 

In the case (-B+) we get a bouncing behavior of scale factor a(r) analogous 
to that in the case (A + ). In the case (Bo) the solution is non-singular, the 
function a(r) is monotonically increasing to infinity from a some non-zero 
value (as in the case (A )). 
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For the (_B + )-case we get 

{\ a C a 2 f <C= (\C 2 ) 2 - K(C 2 )\ (6.15) 

and, hence, (C2) 2 = habC^C^ < 0. The matrix (hob) is not positive-definite. 

Thus, we have shown that the solutions with the bouncing behavior may 
take place only when there are "phantom" fields among the scalar ones. 

For non-singular solutions of the cases (A ) and (B ) we get by analogous 
consideration that (C2) 2 = h a bCS;C\ = 0. In these cases the non-singular 
solutions may take place also for positive-definite matrix (h ab ) and but for 
zero constants: C 2 = 0. 

7 Conclusions 

In this article we have studied a (n + l)-dimensional gravitational model with 
a set of I scalar fields and an antisymmetric form field of rank p + 2. The 
special case of this model containing one scalar field (/ = 1) with positive- 
definite kinetic term was investigated in [T6*j . 

The metric was taken to be diagonal one, and for the (p+ 2)-form a com- 
posite electric Sp-brane form was adopted. All ansatz functions depended 
only on the one distinguished coordinate, u. Under these conditions the ini- 
tial model was reduced to an effective 1-dimensional cx-model. The diagonal 
form of our metric ansatz resulted in the appearance of constraint equations 
on the charge densities of branes associated with the (p + 2)-form field. 

It was found in [TB] that these constraint equations may be satisfied for 
certain odd values of the space-time dimension, given by D = Am + 1 = 
5,9,13..., the model allowed the maximal number of the charged electric 
branes. 

For these special odd dimensions and non-exceptional dilatonic coupling 
(A 2 7^ 4 ( n "_ 1 ) ) we have found exact cosmological-type solutions given in equa- 
tions ()5.H()j) - (J5.32|) . We have singled out a special solution given by ()5.21|) . 
describing an attractor power- law expansion ( see, ()5.H8j) . ()5.38j) and ()5.41j) ). 
e.g. with accelerated expansion of the scale factor. 

We have also found an asymptotical Kasner type behavior of the solutions 
for small (r — > +0) and large (r — > +00) values of proper time when the 
matrix {h^) is positive-definite. In (15.57)) we have presented the relation 
between Kasner parameters (for scalar fields) in both asymptotics. 
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Here we also have found non-singular solutions, e.g. with bouncing be- 
havior of the scale factors. We have shown that solutions with bouncing 
behaviour take place when there are "phantom" scalar fields in the model. 
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